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Abstract 

We examine numerically the distribution function fn(r) of distance r between opposite 
polygonal nodes for random polygons of N nodes with a fixed knot type K. Here we consider 
three knots such as 0, 3i and 3ijj3j.. In a wide range of r, the shape of /ftr(r) is well fitted by 
the scaling form [1] of self-avoiding walks. The fit yields the Gaussian exponents vk = \ and 
jK = 1. Furthermore, if we re-scale the intersegment distance r by the average size Rk of 
random polygons of knot K, the distribution function of the variable t/Rk should become 
the same Gaussian distribution for any large value of N and any knot K. We also introduce 
a fitting formula to the distribution gx{R) of gyration radius R for random polygons under 
some topological constraint K. 

1 Introduction 

Polymer chains in solutions or gels may be highly self-entangled: such entanglements should be 
important to understand some features of polymeric materials. A variety of knots can appear 
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by connecting the two ends of a polymer chain. In fact, various knotted DNAs are synthesized 
in experiments through random closure of nicked DNA chains [2,3]. Since topological ques- 
tions were addressed by Delbriick, Frisch and Wasserman [4,5], several aspects of knotted ring 
polymers, such as the probability of random knotting [6-15], the average sizes [16,17] and the 
complexity of their conformations [18] have been studied numerically and analytically. 

Let us discuss the average size of knotted ring polymers with no excluded volume, i.e., the 
mean-squared gyration radius R 2 K {N) of iV-noded random polygons with fixed knot type K 
[19-25]. We consider random polygons as a simple model of ring polymers in solution at the 9- 
point [26]. At the #-point polymers should have no effect of excluded volume. Furthermore, ring 
polymers keep their topology unchanged. It has now been established in simulations [20,22-24] 
that the average size of random polygons with a fixed knot is larger than that with no topological 
constraint, when N is large. The topological swelling of random polygons may be explained in 
terms of entropic repulsion caused by the topological constraint. The phenomenon should be 
closely related to the "topological excluded volume" proposed for such random polygons that 
possess the trivial knot [19]. Concerning the large- TV behavior of R 2 K {N), however, the 
numerical studies do not unanimously arrive at the same conclusion. We have analyzed the data 
of R 2 K (N) for a model of random polygons [23], assuming the scaling formula of the following 
form: 

R 2 K {N) = A K N 2VK (1 + B k N- Ak + ...) . (1) 

The result favors to the interpretation of vk = fsAW- However, limiting the analysis to a 
narrower range of N, the alternative interpretation vk = ^rw is also possible. Here, self- 
avoiding walks (SAW) and random walks (RW) have the scaling exponent fsAW = 0.588 and 
^RW = 0.5, respectively. Thus, in order to clarify the large-TV behavior of random polygons 
with fixed knots, it would be interesting to investigate some other quantity associated with the 
asymptotic behavior. 

In this paper we study the following two quantities of random polygons with a fixed knot: 
i) the distribution function of the distance between opposite nodes and ii) the distribution of 
the radius of gyration. If the "topological excluded volume" corresponds to a certain amount 
of excluded volume, then the distance between opposite nodes should follow a non-Gaussian 
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distribution. For ring polymers with excluded volume, the distribution of the distance between 
opposite nodes should be non-Gaussian, while it is Gaussian for random polygons. Here we 
assume that the distance between opposite polygonal nodes plays the similar role as the end-to- 
end distance of a linear chain. For the self-avoiding walk, the end-to-end distance distribution 
is non-Gaussian [27,28]. 

Through computer simulation of random polygons with fixed knots, we have evaluated the 
distributions of the distance r between opposite nodes of random polygons under the topological 
constraints [29]. We are concerned with the trivial knot, the trefoil knot, and the composite 
knot consisting of two trefoil knots, which are denoted by 0, 3i, and 3i(J3i, respectively. We 
show that the scaling form of the self-avoiding walk gives good fitting curves to the data of the 
distributions in a wide range. It should be remarkable, since the scaling form of the end-to- 
end distance distribution is derived when p = r/N u is finite and very large [1,27,28,30-34]. 
Furthermore, we show that the distribution function of normalized distance v/Rk should be 
given by the same Gaussian form for any N and K. Here we recall that Rk denotes the average 
size of random polygons of knot K. Thus, it is suggested that the effect of the "topological 
excluded volume" should be different from the standard excluded- volume effect. 

We have also evaluated the distribution of the gyration radius for random polygons under 
some topological constraints [29]. We introduce a formula for describing the distribution, and 
discuss its fitting curves. The formula of the gyration-radius distribution is new, in particular, 
for random polygons under topological constraints. We note that for the Gaussian random walk, 
several approximate formulas of the gyration-radius distribution are known [35,36]. 

The paper consists of the following. In §2, we explain the model of random polygons, and 
define some symbols for the distribution functions. In §3, we describe briefly some procedures 
of the computer simulation. In §4, we discuss the numerical results of the present research. We 
plot the distribution function of the intersegment distance for random polygons under topo- 
logical constraints. Through fitting curves to the data, we discuss that the distribution of the 
intersegment distance should be well approximated by the Gaussian distribution. We also plot 
the distribution of the gyration radius for random polygons under topological constraints. We 
thus investigate topological effects on the average sizes of random polygons. In §5, we discuss 
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that the "topological excluded volume" should be different from the standard excluded volume, 
and present an open question. 

2 Model and distribution functions 

We consider a model of random polygons in which a polygon Vn consists of N line segments 
of length a. It is specified by position vectors of its nodes, Vn = {fi,f2, ■ ■ ■ ,i"n). All cyclic 
permutations of the set of position vectors correspond to the same polygon. We recall that 
random polygons have no excluded volume. Hereafter we set a = 1. 

When a polygon is topologically equivalent to a knot K, we call it a polygon of knot type 
K. The configuration space C of polygons is divided into subspaces Ck in which all polygons 
have the same knot K. We have C = J2k ^k- 

For a polygon Vn, we denote the intersegment vector from the ith node to the i + AiVth 
node 

r(i;X,V N ) =r i+XN -r i7 (2) 

where the progress parameter A takes a value between and 1. Here we assume the convention: 
rN+i = Ti- 

We define the distribution of the distance between the ith and the i + AiVth nodes by the 
probability f(r; A, N)Ar that the length of the intersegment vector r(i;\,VN) takes a value 
between r and r + Ar: 

M N r+Ar 

f(r;\,N)Ar = — Y,Y, / dr 6 (r - \r(i; \,V N , m )\) . (3) 

m=l i=l Jr 

Here Ar is a small positive real number. We choose it so that the statistical fluctuation of 
f(r;X,N) becomes moderately small. The distribution of the distance between two nodes for 
random polygons with a fixed knot type K is similarly defined by 

1 M N r+Ar 

Mr;\,N)Ar = 1 ^ r ^2Yl / dr 5 (r - \r(i; \,V N , m )\) x(V N , m , K), (4) 

m=l i=l r 

Here the indicator function x('P, K) filters the polygons of knot type K; it takes the value 1 if 
V G Ck and otherwise. 
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We calculate the distribution f(r;X,N) of the intersegment distance r by randomly gener- 
ating a large number of polygons VN,m with length TV for m = 1, . . . , M. Here the subscript m 
denotes the mth polygon generated. The number of generated polygons of the knot type K is 
given by M K = £ m x(?V,m, *0, and we have M = £ x M K . 

Let us denote the square of the gyration radius of a polygon Vn by 

R 2 G (V N ) = ^J2(r l -r J ) 2 . (5) 
We define the distribution g(R; N) for gyration radius R by 

i M rR+AR / , \ 

<?(i?; iV) Ai? = tt E / dR5(R - yj ' R 2 G {V N , m ) , (6) 

m=l ■'•R \ / 

and the one for polygons with knot type JT by 

i M pR+AR / , \ 

g K (R;N)AR = — ^ J r dR5[R- 1 R 2 G {V N , m )J x(V N , m ,K). (7) 
3 Simulation procedure 

A pivot move for a polygon is a rotation of a chain of segments, randomly chosen from the 
polygon, around the axis passing the two endmost nodes of the chain by a random amount of 
angle cf) [23,37]. The rotation angle <p is selected randomly from the interval between and 360 
degrees. We do not check self-intersections during the process of rotation of the chain since such 
configurations are negligible in the space C. 

We start from a seed conformation placed on the cubic lattice, which is regarded as a special 
conformation of the off-lattice polygon in the continuum space [38]. We then generate a sequence 
of polygons by applying the pivot moves repeatedly. After discarding the initial 2000 transient 
conformations, we take samples of polygons at every 200 pivot moves. 

To determine the topology of polygons, we employ two simple knot invariants. We calculate 
the special value of the Alexander polynomial Ax{t) at t = — 1 [6] (which is also called the 
determinant of a knot), and the Vassiliev invariant of the second order vi(K) [39,40]. With 
these invariants, the chance of misidentification of topology class for a given polygon should be 
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negligible and much smaller than the statistical errors of the data, as far as the simple knots are 
concerned. 

The simulation has been performed for polygons with the length N = 300 and 600. We have 
generated M = 3 x 10 6 random polygons for each given length N. We have classified those 
polygons into four groups according to their knot types, the three groups of polygons with the 
specific knot types 0, 3i, and 3i(J3i, and the other group of knot types other than the previous 
three. The three knots have distinct sets of values for the two knot invariants |A.r-(£ = — 1)| and 



The distribution function /^(r;A, N) of the intersegment distance r has been evaluated 
at the progress parameter A = 1/4, 1/2 and 3/4 for random polygons under some topological 
constraint K [29]. However, we focus on the case of A = 1/2. The range of intersegment distance 
r is divided into a number of bins of with the width Ar. Here we set Ar = 0.25. We enumerate 
the number of intersegment distances in each of the bins. The distribution function is obtained 
by dividing the number of each bin by the total number of intersegment distances. Similarly we 
numerically evaluate the distribution gxiR] N) of gyration radius R for random polygons under 
some topological constraint K. Here we take AR = 0.25. 

4 Results of the simulation 

4.1 Functional forms of the distributions 

The asymptotic scaling form of the end-to-end distance distribution of the self-avoiding walks 
is derived for the region p = r/N u » 1 [1,30,31]. We now apply it to the data of the distribu- 
tion function for the distance between opposites nodes for random polygons under topological 
constraint K. We thus have the following: 



v 2 {K). 



f K (r;X,N) = 



0k = 



A K r 2+dK exp -D K r 

du K + 1 - j K - d/2 
l-u K 



(8) 



(9) 



8 K = 



\-v K 



(10) 



Hereafter we set d = 3. 
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For the distribution gx (R; N) of gyration radius R, we introduce the following formula: 

gK{R;N) = Ag >K \R-c K \ e *Kexp[-Dg,K\R-c K \ s » K ] ■ (11) 

For the Gaussian random walk there are some approximate expressions for the distribution of 
the gyration radius [35,41]. (See also §8 of Ref. [36].) For instance, the large R case of Fixman's 
result [35] corresponds to a special case of the formula where we set 5 9: k = 2, ck = 0, and 

0g,K = 1- 

4.2 Distribution function fx{r; X, N) of intersegment distance r 

The intersegment distributions fx^r; A, N) at A = 1/2 for N = 300 and 600 are presented 
in figures Q and [21 respectively. Here the topological conditions denoted by K correspond to 
restriction of random polygons into the following sets: (i) all polygons; (ii) polygons of the 
trivial knot 0; (iii) polygons of the trefoil knot 3i; (iv) polygons of the composite knot 3i(J3i; 
(v) polygons of any knot types other than the three knots 0, 3i, and 3ijJ3i. We thus consider 
the five different topological conditions. We note that the case (i) corresponds to no topological 
constraint. We denote the distribution functions of the five cases simply as f a u, fy, f^, f3 1 ^3 1 , 
and fathers, respectively. 

The fitting curves of figures ^ and [2] are fit well to the data points. The curves are given by 
the scaling form ((HI), and are all very close to the Gaussian distributions. Here we note that it 
is also the case with the data for A = 1/4 and 3/4. The numerical estimates for the exponents 
9k and 8 k and the constants Ak and Dk are given in table ^ The actual ranges of distance r 
used for the fitting curves are also shown in tabled The fitting curves fit very well to the data 
points not only in the range of r larger than the peak position but almost in the entire range of 
r. The x 2 values per datum are very small. Very small deviations are only seen in the small r 
region, although the region is out of the fitting ranges. 

The best estimates of the exponents 9k and 5k almost agree with the Gaussian values, i.e., 
vk ~ 1/2 and jk ~ 1 5 within the range of estimation errors, for all the five different topological 
conditions and for both N = 300 and 600. The constant Dk depends on the polygonal length 
N. However, it does not change very much for the different knot types with respect to the 



7 



estimation errors. The constant Ak depends on the knot type K for N = 300. However, the 
difference among Ak's becomes smaller for N = 600 than for N = 300. It is thus suggested 
that they should become the same value when N is very large. 

Let us denote by r* K the peak position of the distribution /^(r). Assuming the scaling form 
(jHJ, the peak position r* K is given by 

The peak position r* K may characterize the knot dependence of the distribution function fxif)- 
When uk = 0.5 and jk = 1-0, the form of fK( r ) is determined by the parameter Dk- 

In figure^ the peak position of the distribution f$(r) is larger than r* u of f a ii- For fothers, 
r others ^ s smaller than r* all . In figure El the peak positions of f^, and f3 1 ^ 1 are all larger 
than that of f a u for N = 600. Their values of N = 600 are much closer to each other than in 
the case of N — 300. Here the peak position ^Q^ ers of /others is smaller than that of f a n also in 
the case of N = 600. It is thus suggested that when N is very large, the peak positions of //<-(r) 
of simple knots should be given by the same value and the distributions fxir) should approach 
a universal form. 

The observations in figures^and EJsuggest that fixing a knot type of a random polygon leads 
to effective repulsion or attraction among internal segments of the polygon depending on the 
complexity of the knot type. When the length iV becomes very large, polygons of very complex 
knots can appear. They should have smaller conformations than other polygons of simpler 
knots. As we see in figure ^for N = 300, random polygons with the trivial knot have larger 
conformations on the average than those of no topological constraints, while random polygons 
of more complex knots have smaller conformations. This should be consistent with the effective 
swelling observed in the studies on the average sizes of random polygons with some fixed knots 
[20,22-24]. 

Let us discuss the A- and iV-dependence of the distribution function fx(r,\,N) for a knot 
K. We denote by rx(A,A r ) the average of the intersegment distance r at the parameter A for 
random polygons of iV nodes with the knot K. We shall suggest that for a given knot K, the 
distribution function /^(r, A, N) should depend on N and A only through the value rx(A, N). 
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We introduce the distribution fx of normalized intersegment distance x = r/rx- Here we 
note fx(x, A, N)dx = fx(f, X, N)dr. The data for the three knots show that the function fx 
does not depend on either A or N. In figure |3] the data points of the distribution function 
f$)(x; X,N) of the normalized intersegment distance x = r/r^ are shown for the four cases: A = 
1/2 or 1/4 and N = 300 or 600. It is clear that the data points for all the four cases are located 
on the same curve. The best estimates of the fitting parameters to the data of fx are given in 
tables El and |31 for A = 1/2 and 1/4, respectively. As far as the five topological conditions are 
concerned, each of the fitting parameters of a condition K has almost the same value for A^=300 
or 600 and for A = 1/2 or 1/4. Thus, we suggest that for a given knot K the distribution of the 
normalized intersegment distance, fx(x,X,N), should be given by the same Gaussian form for 
any N and A. 

We now show that the A- and iV-dependence of the average distance rx(X, N) is given by the 
Gaussian one in the cases of A = 1/4, 1/2 and 3/4 for the three knots. Let us denote by P(r; N) 
the end-to-end distance distribution of the Gaussian random walk of N steps. For random 
polygons consisting of two Gaussian chains of AiV steps and (1 — X)N steps, the probability 
distribution of the vector r connecting the two end-points should be proportional to the product 
-P(r; XN)P(r; (1 — X)N). Here we note that the intersegment vector r for the parameter A 
connects the common end-points of the two random walks of length AA?" and (1 — X)N. For the 
case of no topological constraint, therefore, the constant D a ii is given by 

3 

DaLl = 2A(1 - A)iVa 2 (13) 
Here a denotes the bond length, and a = 1 in the paper. We thus have 



r all (X,N) = J—^X(l-X)Na (14) 

V o7T 

The ratio rx(X,N)/r a u(X,N) is plotted in figure |I] for N = 300 and 600 with respect to the 
topological conditions, 0, 3i, 3ijJ3i, and others. In each case of the four topological conditions, 
the ratio takes almost the same value for A = 1/2, 1/4 and 3/4. Furthermore, we find that 
the ratio should also coincide with the ratio of the gyration radii, Rx/R a ii- Thus, we have the 
conjecture: rx(X,N) = r a u(X, N)Rx / Rail- If the conjecture is true, then the A-dependence of 
the distribution fx(r, X, N) is completely given by the Gaussian one, and the A^-dependence is 



given by the Gaussian with the rescaling factor Rk / R a ii ■ Here we remark that the ratio Rk / Rail 
may depend on the number of nodes N, since random polygons under a topological constraint 
can be larger or smaller than that of no topological constraint due to entropic repulsion induced 
by the topological constraint. 

We explain some details about the average size of random polygons. We first recall that the 
symbol Rk denotes the square root of the mean square radius of gyration for random polygons 
of a topological constraint K. We denote by (Rk) the mean gyration radius of a polygon 
averaged over an ensemble of random polygons of a topological constraint K. In figure 01 the 
ratio (Rk) / {Rail) is plotted for the four topological conditions. However, the difference between 
the two ratios, Rk/RuII and (Rk) / '{Rail) > should be smaller than the error bars. 

Let us discuss the knot-dependence of the distribution fx- We show that the distribution 
fx of normalized distance x for a knot K should be almost independent of the knot type. In 
figure E] the rescaled distribution fx for the five topological conditions are plotted against the 
normalized intersegment distance x = r/rx for the case of N = 300 and A = 1/2. We see in 
figureEJthat the distributions fx should almost the same for the three knots. It is consistent with 
the observation that the estimates of the parameters are of similar values for the five topological 
conditions, as shown in tables El and El Thus, the knot-dependence of the distribution fx should 
be renormalized into the value r^(A, N). 

4.3 Distribution g K (R;N) of gyration radius R 

The distribution functions gx(R',N) of gyration radius R are shown in figures E] and [7| with 
respect to the five topological conditions for N = 300 and 600, respectively. Here the five cases 
are the same as in §4.2. We denote the distributions gx for the five cases briefly as g a a, g$, 
93d 931(131) an d g thersi respectively. The fitting curves in figures El and fit well to the data 
points in some ranges of R. The curves are given by the formula (|11|) . Thus, the formula 1)11(1 
approximates the distribution gx(R', N) of gyration radius effectively. They could be useful for 
studying topological effects on the gyration radius. 

Let us consider plotting the distribution gx with respect to a normalized variable of the 
gyration radius, R/Rx- More precisely, we consider plotting gx in terms of the variable Rj (Rx)- 
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Here we recall that {Rk) denotes the average of the gyration radius of a polygon averaged over 
an ensemble of random polygons with a given knot K. We introduce distribution gK(y,N) of 
the variable y = R/{Rk) by the relation: g^iy; N)dy = gx(R', N)dR. Hereafter, however, we 
denote (Rr) simply by Rk, except for figure captions. 

We now present the rescaled distribution gx in figures |H1 and E3 for N = 300 and N = 600, 
respectively. We find that for a given knot K the distribution gx of normalized gyration radius 
y = R/Rk should be independent of the knot type almost completely. In figures El and the 
data points and their fitting curves of gx overlap each other for the cases of the three knots, 0, 
3i, and 3ijJ3i. We also find that the fitting curves are fit well to the data points both in figures|H] 
and|§J They are drawn by a fitting formula corresponding to (1111) 

(15) 

The fitting parameters A g ,K, D g K and ck correspond to A 9j k, ^g,K and ck of the formula (fTTj) 

as 

A 9: k = Ag^K/R 1 ^ 9 '* j D g,K = Dg tK /R S £ K , c K = c K Rk ■ (16) 

The best estimates of the fitting parameters are listed in table @] The x 2 values are good, in 
particular, for the cases of the three knots, g^, g^ and ■ Here the fitting range of y = R/Rk 
is given by from 0.4 to 2.0 for all the fitting parameters given in tabic 21 

The knot-dependence and the iV-dependence of distribution gK should be renormalized into 
the mean square radius of gyration, R K (N). The rescaled distributions c/k for the three knots 
do not depend on the polygonal length, N. The fitting curves of c/k for the three knots are 
almost the same for N = 300 and N = 600. We can confirm the observation in figures |H1 
and|2lby comparing the estimates shown in table 01 The fitting parameters 0g t K, $g,K, A 9t K 
and ck depend on neither the knot type K nor the polygonal length N with respect to their 
errors. The normalization of gyration radius, R/Rk, should be thus essential when analyzing 
the distribution function gK- 

Random polygons of relatively simple knots should be larger in size than the average one, 
while those of more complex knots should be smaller. It depends on the polygonal length N 
whether the size of random polygons of a given knot should be larger or smaller than the average. 
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gxiy; N) = Ag jK \y - c K \ e ^ K exp 



-D. 



g,K\y - C K\ 



In figure|13 the peak position of the distribution is larger than those of the other distributions 
9all, 531, 53x831 and gathers- In figure the peaks of g 9 , g 3l , g 3li3l are clearly located on the 
right hand side of the peak of g a u, while the peak of g thers is located on the left hand side. The 
equilibrium length of a random knot can be a criterion whether it is larger or smaller than the 
average [42] . We also note that the peak positions of distributions gx shown in figures El and 
are roughly consistent with the average values of Rx [20,22-24]. 

5 Discussion 

We have found that the distribution fx of intersegment distance for random polygons under 
topological constraint K is almost given by the Gaussian distribution. Furthermore, rescaling 
the distance by the average distance rx, we have shown that the A- and ^-dependence of fx 
is renormalized into the average distance nr(A, N). We have proposed the conjecture: rx = 
r aii Rfc/Raii, for any A, N and K. Here we have assumed that ./V is large enough. If it is true, 
then topological constraints do not have any effect on the distribution of intersegment distance, 
fx, except for scaling the distance by the factor Rx/R a u- 

The effect of the "topological excluded volume" should be rather different from the standard 
excluded volume effect of self- avoiding walks. It does not correspond to a real excluded volume, 
although the ratio Rx/R a ii of a knot K can become larger than 1 in the case of large N. 
When N increases, random polygons with more complex knots can appear, which should be 
smaller than those of a simple knot. If we consider only such random polygons that have a fixed 
simple knot, then the size can be larger than the average one when N is very large. Topological 
constraints thus may induce effective swelling of random polygons. However, they do not change 
the functional form of the distribution of the distance between two segments. 

Let us discuss the difference in terms of critical exponents. We denote by v' K the scaling 
exponent defined for the asymptotic behavior of the average size of SAW such as given in fT]). 
For SAW, the exponent v' K corresponds to the exponent vx determined by the formula (JHJ) [1]. 
If des Cloizeaux's relations © and (jlOj) could be valid for random polygons under topological 
constraints, we should have v' K ~ 0.50 from the best estimates for the distribution fx(j", A,iV) 
as shown in table 
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Within the scope of the present research, however, it is not clear whether the two exponents 
vk and v' K should be equal or not. Moreover, it is not clear whether the relations Q and 
(|10jl should be valid for random polygons under topological constraints. It seems that the form 
of the distribution fjc(r;X,N) remains Gaussian with the exponent i/# ~ 0.50 in the limit 
N — > oo. However, the average size R? K (N) might follow the scaling form with a different 
exponent, v' K > 0.5. 
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Figure 1: Distribution /ic(r; A, N) of intersegment distance r at A = 1/2 for N = 300. Dots(.), 
crosses(x), double crosses(*), open squares(D) and closed squares(B) denote the plots of condi- 
tions, all, 0, 3i, 3it)3i, and others, respectively. The plots and fitting curves for all, 0, 3i, 3ijJ3i, 
and others are colored with red, blue, fuchsia, black and orange, respectively. 
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Figure 2: Distributions ficfr) \ N) of intersegment distance r at A = 1/2 for N = 600. Dots(.), 
crosses(x), double crosses(*), open squares(D) and closed squares(H) denote the plots of con- 
ditions, all, 0, 3i, 3ijJ3i, and others, respectively. The plots and fitting curves are colored with 
red, blue, fuchsia, black and orange, respectively. 
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Figure 3: Distributions f$(x; A, iV) of normalized intersegment distance x = r/rx at A = 1/2, 1/4 
and for N = 300,600. Dots(.), crosses(x), double crosses(*) and open squares(D) denote the 
plots with (A, iV) = (1/2, 300), (1/4, 300), (1/2, 600) and (1/4,600), respectively. They are dis- 
played with red, blue, fuchsia and black, respectively. Here Ar of eq. Q is given by Ax = 0.02. 
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Figure 4: The ratio of the average distances rR/r a u and that of the gyration radii (Rk) / (Rail) 
for a given topological condition K. The ratios vr(X, N)/r a n(X, N) for A = 1/2, 1/4, and 3/4 
are denoted by dots(.), crosses(x), and double crosses(*), respectively. Here K is given by 0, 
3i, 3ijj3i, and others {oth). The ratio (Rr) I '(Rail) is denoted by open squares(D) for the four 
cases of K. Here we consider N = 300 and 600. Dots, crosses, double crosses and squares are 
colored with red, blue, fuchsia and black, respectively. 
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Figure 5: Distribution fx(x; \,N) of normalized intersegment distance x = r/rx for A = 1/2 
and N = 300. The distributions fK(x;\,N) for all, 0, 3i, 3iJJ3i, and others are denoted by 
dots(.), crosses(x), double crosses(*), open squares(D) and closed squares(H), respectively. The 
plots and fitting curves are colored with red, blue, fuchsia, black and orange, respectively. Here 
Ar of eq. @ is given by Ax = 0.02. 
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iV=300 




Figure 6: Distribution gx{R;N) of gyration radius R for N = 300. Dots(.), crosses(x), double 
crosses(*), open squares(D) and closed squares(H) denote the distribution gx(R;N) for all, 
0, 3i, 3i {}3i and others, respectively. The plots and fitting curves are colored with red, blue, 
fuchsia, black and orange, respectively. 
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Figure 7: Distribution gx{R;N) of gyration radius R for N = 600. Dots(.), crosses(x), double 
crosses(*), open squares(D) and closed squares(H) denote the distribution gK{R;N) for all, 0, 
3i, 3i)J3i, and others, respectively. The plots and fitting curves are colored with red, blue, 
fuchsia, black and orange, respectively. 
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Figure 8: Distribution gx{y]N) of the normalized gyration radius y = R/(Rk) for N = 300. 
Dots(.), crosses(x), double crosses(*), open squares(D) and closed squares(H) denote the dis- 
tribution gx{y',N) for the conditions, all, 0, 3i, 3ij)3i, and others, respectively. The plots and 
fitting curves are colored with red, blue, fuchsia, black and orange, respectively. Here AR of eq. 
(0 is given by Ay = 0.02. 
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iV=600 




Figure 9: Distribution g^-(y;iV) of the normalized gyration radius y = R/Rk for N = 600. 
Dots(.), crosses(x), double crosses(*), open squares(D) and closed squares(H) denote the dis- 
tribution gjc(y;N) for the conditions, all, 0, 3i, 3ijJ3i, and others, respectively. The plots and 
fitting curves are colored with red, blue, fuchsia, black and orange, respectively. Here AR of eq. 
(0) is given by Ay = 0.02. 
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Table 1: The fitting values of the scaling formula (jHJ) to the data of the distribution /^(r; A, N) 
at A = 1/2 with \ 2 values per datum. 



N 


K 


IK 


vk 


D K x 10 2 


A K x 10 3 


x 2 


Fitting 
Range 


300 


all 


1.04 ±0.03 


0.503 ±0.005 


1.9 ± 0.1 


1.78 ±0.09 


1.13 


5.75 - 20 







0.87 ±0.13 


0.509 ±0.019 


1.7 ±0.5 


0.8 ±0.1 


1.05 


6.5 - 16 




3i 


0.93 ±0.09 


0.50 ±0.01 


1.9 ± 0.3 


1.1 ±0.1 


0.90 


6.25 - 19 




3i#3i 


1.0 ± 0.1 


A r 1 I A AO 

0.51 ± 0.02 


OA I A C 

2.0 ± O.o 


1.5 ± 0.3 


A O K 

0.85 


5.75 — 17 




others 


1.1 ± 0.1 


0.48 ±0.02 


2.9 ±0.6 


4.0 ±0.5 


0.81 


5.5-15 


600 


all 


1.03 ±0.02 


0.501 ±0.004 


0.99 ±0.05 


0.63 ±0.03 


0.76 


7.5 - 30 







0.8 ±0.2 


0.51 ±0.03 


0.8 ±0.4 


0.19 ±0.07 


0.78 


9-23 




3i 


0.8 ±0.2 


0.51 ±0.02 


0.9 ±0.3 


0.22 ±0.06 


1.22 


8.5-23.75 




3i#3i 


1.0 ±0.1 


0.52 ±0.02 


0.7 ±0.3 


0.29 ±0.06 


0.92 


7.5-22 




others 


1.06 ±0.05 


0.493 ±0.008 


1.2 ±0.1 


0.94 ±0.09 


0.97 


8.25 - 28.5 



Table 2: The fitting values of the scaling formula (jHJ) to the data of the re-scaled distribution 
A, N) at A = 1/2 with x 2 values per datum. 



Fitting 



N 


K 


IK 


VK 


D K 


A K x 10 2 


x 2 


Ran 


ge 


300 


all 


1.04 ±0.04 


0.504 ±0.006 


1.25 ±0.04 


6.3 ±0.2 


0.80 


0.70- 


2.00 







0.89 ±0.08 


0.51 ±0.01 


1.39 ±0.08 


7.7 ±0.6 


1.20 


0.70- 


2.00 




3i 


0.99 ±0.08 


0.52 ±0.01 


1.30 ±0.08 


6.9 ±0.5 


0.89 


0.70- 


2.00 




3i#3i 


1.1 ±0.1 


0.52 ±0.02 


1.2 ±0.1 


6.4 ±0.8 


0.90 


0.70- 


2.00 




others 


1.04 ±0.07 


0.49 ±0.01 


1.26 ±0.07 


6.2 ±0.4 


0.91 


0.70- 


2.00 


600 


all 


1.03 ±0.04 


0.502 ±0.006 


1.26 ±0.04 


6.4 ±0.2 


0.75 


0.70- 


2.00 







1.0 ±0.1 


0.53 ±0.02 


1.3 ±0.1 


7.2 ±0.8 


0.81 


0.70- 


2.00 




3i 


0.9 ±0.1 


0.51 ±0.02 


1.4 ±0.1 


7.8 ± 1.0 


1.26 


0.70- 


2.00 




3i#3i 


1.0 ±0.1 


0.52 ±0.02 


1.3 ±0.1 


7.2 ± 1.0 


1.05 


0.70- 


2.00 




others 


1.06 ±0.05 


0.499 ±0.008 


1.24 ±0.05 


6.2 ±0.3 


1.00 


0.70- 


2.00 
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Table 3: The fitting values of the scaling formula (JBJ to the data of the re-scaled distribution 
fx(x; X, N) at A = 1/4 with x 2 values per datum. 



Fitting 



N 


K 


IK 


VK 


D K 


A K x 10 2 


x 2 


Ran 


ge 


300 


all 


1.10 ±0.04 


0.514 ±0.006 


1.19 ±0.04 


6.0 ±0.2 


0.88 


0.70- 


2.00 







1.01 ±0.09 


0.52 ±0.01 


1.28 ±0.08 


6.8 ±0.6 


1.54 


0.70- 


2.00 




3i 


1.15 ±0.07 


0.54 ±0.01 


1.14 ±0.07 


5.9 ±0.4 


0.94 


0.70- 


2.00 




3i#3i 


1.0 ±0.2 


0.51 ±0.02 


1.3 ±0.2 


6.8 ± 1.0 


1.06 


0.70- 


2.00 




others 


1.11 ±0.07 


0.50 ±0.01 


1.20 ±0.07 


5.8 ±0.4 


1.03 


0.70- 


2.00 


600 


all 


1.05 ±0.05 


0.507 ±0.007 


1.24 ±0.05 


6.2 ±0.3 


1.28 


0.70- 


2.00 







1.0 ±0.1 


0.53 ±0.02 


1.3 ±0.1 


7.0 ±0.7 


0.78 


0.70- 


2.00 




3i 


0.8 ±0.1 


0.50 ±0.02 


1.4 ±0.1 


8.1 ± 1.1 


1.24 


0.70- 


2.00 




3i#3i 


0.9 ±0.2 


0.51 ±0.02 


1.4 ±0.1 


7.5 ± 1.1 


1.18 


0.70- 


2.00 




others 


1.13 ±0.05 


0.511 ±0.008 


1.17 ±0.05 


5.8 ±0.3 


1.07 


0.70- 


2.00 



Table 4: The fitting values of the formula (|15|) to the data of the re-scaled distribution qk{v\ N) 
of the normalized gyration radius y = R/(Rk) for iV-noded random polygons of topological 
condition K with x 2 values per datum. The fitting range of y is from 0.4 to 2.0. 



N 


K 






D g ,K 


A giK x 10- 4 


CK 


x 2 


300 


all 


6.2 ±0.3 


1.44 ±0.03 


11.0 ±0.4 


0.019 ±0.009 


0.423 ±0.007 


3.63 







7.9 ±0.4 


1.31 ±0.03 


14.6 ±0.6 


0.4 ±0.3 


0.441 ±0.006 


1.61 




3 X 


8.3 ±0.3 


1.21 ±0.03 


16.1 ±0.6 


1.7 ± 1.1 


0.450 ±0.004 


3.33 




3i#3i 


8.7 ±0.5 


1.13 ±0.03 


17.3 ±0.8 


6±5 


0.457 ±0.005 


1.16 




others 


7.7 ±0.3 


1.12 ±0.03 


15.1 ±0.6 


1.1 ±0.6 


0.441 ± 0.004 


3.97 


600 


all 


6.2 ±0.3 


1.42 ±0.02 


11.0 ±0.3 


0.020 ±0.007 


0.424 ±0.006 


5.91 







7.7 ±0.3 


1.30 ±0.03 


14.9 ±0.4 


0.46 ±0.23 


0.457 ±0.005 


1.86 




3i 


8.3 ±0.3 


1.21 ±0.03 


16.3 ±0.5 


1.9 ±1.1 


0.458 ± 0.004 


2.39 




3i#3i 


8.1 ±0.3 


1.22 ±0.03 


16.0 ±0.5 


1.3 ±0.6 


0.463 ± 0.003 


2.38 




others 


7.0 ±0.3 


1.25 ±0.03 


13.1 ±0.5 


0.14 ±0.08 


0.434 ±0.006 


6.90 
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